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ABSTRACT 


The  seismic  data  analyzed  include  an  Aleutian  earthquake, 
several  samples  of  ambient  noise  along  with  the  noise  occur- 
ing  before  the  event,  and  a  25  second  sample  of  coda  follow¬ 
ing  the  signal.  The  signal  coda  sampled  had  nearly  uniform 
power.  A  preliminary  interpretation  of  the  noise  and  signals 
strongly  suggests  the  importance  of  energy  conversions  possibly 
due  to  the  complex  geology  at  APOK.  The  geology  is  known  to 
be  very  complex  and  is  characterized  to  first  order  a  thick 
wedge  of  low  velocity  strata  imbedded  in  material  of  higher 
propagation  velocity.  The  evidence  of  conversion  is  based 
on  the  asymmetry  of  the  F-K  spectrum  of  ambient  noise  with 
respect  to  positive  and  negative  wave  number.  The  results 
indicate  a  predominance  of  obliquely  incident  up-going  waves 
at  all  frequencies.  Further,  the  signals  and  coda  similarly 
displayed  anomalous  down-going  pulses  attenuated  by  approx¬ 
imately  6  db  with  respect  to  the  up-going  pulse.  Further, 
the  down-going  pulse  contains  lower  frequencies  than  the 
up-going  pulse  and  has  a  much  lower  apparent  vertical  velocity. 
The  results  obtained  by  mirror  imaging  the  vertical  array 
indicated  that  a  loss  of  approximately  3  db  would  be  encount¬ 
ered  in  using  a  conventional  signal  model  at  this  site. 


ARRAY  RESPONSE 


Due  to  the  finite  dimension  of  the  array  and  the  uneven 
spacing  of  sensors,  and  finite  duration  of  the  samples,  the 
response  to  a  sinusoidal  wave  input  can  be  very  complicated 
and  can  possibly  lead  to  misinterpretations.  For  an  analysis 
of  the  response  of  a  continuous  finite  array  consider  the  re¬ 
sponse  to  sine  wave  inputs. 


p(t,x)  \ 


when  and 


exp[iu>o (t-x/c )  j  +  exp[iu>o(t+x/c)]  X^xfXj 


0 


elsewhere 


X2  T2 


F(to.,k) 


.  -itot  -ikx  , 
p(t,x)  e  e  dtdx 


S (to)  Y(k) 


S  (to ) 


1  r  „i(to„-w)T0  „i(to„-<o)T.1  ] 

l  e  0  2  'e  °  1 J 

o 


PT(to)  =  S(to)  STtoT 


2.  1 1-cos  £  (to^-w )  (T2  -  Ti^]} 

(to  -to)2 
o 


Let  Y (k)  =  Y1  (<i>0/c  -  k)  +  Y2  (toQ  +  k) 


Y,  (k) 


i(to  /c  -  k) 
o 


i(to  /c  -  k) 
e  o 


X„  i(to  /c  +  k) 
2  -e  o 


?  * 
k2 


1 


’  (k)  =  - i -  [ei(“o/c  +  k)  X2  -ei(“o/o  +  k)  Xll 

1  i(a>  /c  +  k)  L  J 


PX(k)  =  (Y1  +  Y2)  (Y1  +  Y2)  =  Y1Y1  +  Y2Y2  +  (Y1Y2  +  Y2Y1) 


Y1Y1 


2.  {l  -  cos  [(a>0/c  -  k)  (X2  -  X1)]} 


( in  /c.  -  k )  2 
o 


2‘  I1  -  COS  [(uo/G  +  k)  (X2  "  V]} 
2  2  - 


C oj  /c  +  k)2 
o 


Y1Y2  *  Y2Y1 


(uo/c  -  k)  (wq/c  +  k) 


-  cos  I (w  / c  -  k)  X,  - 
to  1 


-  cos  (w  /c  -  k) 
L  O 


X.  - 


{cos  2kX2  +  cos  2kX1 

(o>o/c  +  k)  X2] 

(u0/c  +  k)  Xj} 


The  total  power  response  of  the  array  is 
R  (u)  =  PT  (w)  Px  (k) 

From  this  it  is  clear  that  the  response  of  the  F-K  spectrum 
to  a  sinusoidal  wave  of  given  frequency  and  wave  number  is  a 
complicated  function  of  both  frequency  and  wave  number  and  in 
practice  is  even  much  more  so  due  to  finite  and  uneven  samp¬ 
ling  in  space.  Thus  to  make  rigorous  use  of  the  F-K  analysis 
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as  a  tool  one  should  compare  a  model  with  the  observed  data, 
and  conclude  that  the  model  is  adequate  only  if  the  F-K  spec¬ 
tra  qualitatively  match  the  model.  In  this  preliminary  report, 
this  is  not  stressed  enough,  but  as  a  cautionary  aid  the  re¬ 
sponse  of  an  impulse  for  the  case  of  infinite  apparent  vertical 
phase  velocity  is  shown  with  each  F-K  spectrum,  and  several 
synthetic  cases  are  run  for  simple  models  of  a  signal  and  a 
model  of  the  ambient  noise. 

The  signal  model  is  generated  with  1.25  cps  pulse  with 
a  .  8  second  echo  at  the  source  and  with  a  receiver  echo  de¬ 
layed  by  using  appropriate  uph.ole  times  obtained  from  pro¬ 
pagation  velocities  observed  at  APOK  and  with  a  surface  re¬ 
flection  coefficient  of  0.9.  The  signal  model  F-K  spectrum 
is  shown  on  Figure  1.  The  split  peak  in  the  spectrum  is  due 
to  the  source  echo  which  nulls  at  1.25  cps.  The  strip  at  the 
bottom  shows  the  array  response . 

The  noise  is  simulated  by  taking  random  numbers  from  a 
Gaussian  population  and  passing  them  through  a  tuned  filter 
at  .25  cps  and  2.0  cps  to  obtain  a  model  of  the  noise  at  the 
surface.  The  model  of  the  noise  at  underlying  depth  is  ob¬ 
tained  from  a  stationary  Markov  chain;  for  example,  the 
noise  at  the  ith  level  is  taken  as  a  fraction  of  the  noise  of 
the  (i  --  l)^*1  channel  added  to  a  new  random  realization 
passed  through  the  tuned  filters.  An  example  of  this  noise 
model  is  shown  on  Figure  2,  It  is  similar  to  ambient  noise 
observed  at  APOK  (see  Figure  3),  except  that  the  highly  cor¬ 
related  noise  peaks  observed  in  the  signal  band  ( .  1<_ f <2 . ) 
were  not  put  into  the  model.  In  the  model  the  sharp  spec¬ 
tral  peak  at  .25  cps  (Figure  2)  is  for  highly  correlated 
noise  between  channels  contrasted  with  that  at  2.0  cps  where 
the  noise  which  is  uncorrelated  between  channels,  with  the 
results  that  the  peak  is  spread  broaaly  over  all  wave  numbers. 
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AMBIENT  NOISE 


The  ambient  noise  spectrum  derived  from  a  four  minute 
sample  is  shown  on  Figure  3.  Comparing  this  with  noise  gen¬ 
erated  from  a  Markovian  process ,  the  .25  and  2.0  cps  peaks 
are  of  similar  character  suggesting  very  high  correlation 
between  channels  for  the  .25  cps  peak  and  very  low  corre¬ 
lation  for  the  2.0  cps  noise  peak.  The  principal  difference 
between  the  observed  noise  at  APOK  on  Figure  3  and  the  syn¬ 
thetic  noise  generated  using  the  extremely  simple  linear  state 
model  is  a  rotation  of  the  whole  pattern  toward  negative  wave 
numbers.  This  same  effect  can  be  produced  by  inputing  the 
random  function  to  a  process  which  produces  negative  delays 
or  lead  time  equal  to  X/C  representing  conversion  to  up-going 
waves  where  the  apparent  vertical  phase  velocity  C  is  obtain¬ 
ed  from  the  slope  of  the  line  shown  on  Figure  3.  The  value 
obtained  for  C  is  approximately  12  km/sec  corresponding  to  an 
incidence  angle  of  about  75°.  This  suggests  the  possibility 
of  Stonely  waves  guided  upwards  along  the  thick  low-velocity 
layer,  dipping  15°.  This  possibility  is  qualitatively  con¬ 
sistent  with  the  anomalous  signal  shown  on  Figure  5. 

Other  differences  between  observed  noise  on  Figure  3 
and  the  model  on  Figure  2  are  the  three  noise  peaks  at  1.0 
cps,  1.4  cps,  and  1.6  cps.  The  1.0  and  1.4  peaks  appear  to 
be  highly  correlated  be. tween  channels;  the  1.6  shows  low 
correlation  in  the  noise  between  adjacent  channels.  These 
peaks  in  the  signal  band  appear  to  have  nearly  infinite 
vertical  phase  velocity  and  are  probably  due  to  Rayleigh 
waves,  i.e.,  vertical  and  possibly  also  horizontal  standing 
waves  trapped  in  the  basin  bounded  by  higher  velocity  base¬ 
ment  complex  rocks. 


SIGNAL,  NOISE  PRECEDING  SIGNAL,  AND  CODA 


A  30-second  noise  sample  before  the  arrival  of  the  Aleu¬ 
tian  event  is  shown  on  Figure  4.  Comparing  this  with  the  am- 
ient  noise  sampler  on  Figure  3  we  note  the  same  asymmetry 
in  the  F-K  spectrum.  In  addition,  up-going  body  waves  are 
indicated  by  peaks  at  .85  cps,  1.25  cps,  1.5  cps  and  2.0  cps. 
There  is  much  weaker  indication  of  down-going  waves  at  .85 
cps,  down  by  3  db  from  the  up-going  waves.  Waves  of  infinite 
vertical  phase  velocity  (possibly  Rayleigh  waves,)  are  indi¬ 
cated  by  the  series  of  peaks  occurring  along  the  frequency 
axis  at  k  =  0.  For  this  sample,  these  peaks  are  indicated 
with  approximately  the  same  power  as  the  up-going  body  waves. 

A  6-second  sample  of  the  earthquake  pulse  is  shown  on 
Figure  5.  The  up-going  pulse  gives  spectral  peaks  at  .85  cps, 
1.20  cps,  and  1.9  cps.  The  apparent  vertical  phase  velocity 
is  'approximately  the  same  as  that  shown  by  Figure  4  for  the 
noise  preceding  the  signal.  Lower  than  expected  vertical 
phase  velocities  suggest  departure  from  the  simple  model  of 
a  pulse  and  echo  based  on  acoustic  log  velocities  (Figure  1). 
The  apparent  velocities  are  lower  by  at  least  fifteen  to 
twenty  percent.  Also,  the  down-going  earthquake  pulse  is 
even  more  anomalous.  The  amplitude  is  down  6  db  from  that 
of  the  up-going  pulse;  the  apparent  vertical  velocity  is  very 
low;  and  the  .85  cps  peak  down-going  phase  appears  to  contain 

lower  frequency.  A  possible  explanation  ofv'the  anomalous  sig- 

\ 

nal  can  be  based  on  dipping  beds. 

This  may  help  to  explain  the  anomalous  lov:  amplitude 
down-going  reflection.  The  anomalous  apparent  vertical  velo¬ 
cities  may  result  from  forward  scattered  P-S  conversions, 
especially  at  the  surface,  due  to  anomalously  high  angle  of 
emergence.  Looking  again  at  Figure  5,  there  appears  to  be 
signal  peaks  at  nearly  infinite  vertical  phase  velocity. 

-  5  - 


Although  possibly  due  to  the  array  response,  no  such  effect 
is  observed  on  the  simulated  signal  on  Figure  1  which,  of 
course,  uses  the  same  array  geometry. 

A  30-second  sample  of  the  coda  following  the  signal  is 
shown  on  Figure  6.  These  show  apparent  spectral  peaks  in 
the  signal  at  .8  cps,  1.1  cps ,  and  1.35  cps .  The  overall 
character  of  the  coda  F-K  spectrum  is  more  similar  to  the 
signal  than  the  noise  preceding  the  signal,  but  is  yet  con¬ 
siderably  different  in  detail  from  that  of  the  signal. 

Figures  7,  8,  9,  and  10  are  processed  by  taking  the 
mirror  image  of  the  vertical  array,  which  aligns  the  up- 
going  and  down-going  pulses  into  a  single  step-out  pattern. 
This  will  attenuate  conversions  to  up  or  down-goind  P  or  S 
waves  and  amplify  normally  reflected  P  pulses  and  Rayleigh 
waves.  It  will  also  effectively  double  the  aperture  of  the 
array  for  these  kinds  of  waves.  Figure  7  for  the  ambient 
noise  shows  body  waves  down  3  to  6  db  from  those  of  appar¬ 
ently  infinite  vertical  phase  velocity.  Similar  results 
are  shown  for  the  noise  before  the  Aleutian  event  on  Figure 
8,  also  showing  the  noise  field  to  be  down  considerably 
lower  in  the  signal  event.  However,  due  to  the  grossly 
anomalous  nature  of  the  signal  at  APOK,  the  signal  and  coda 
are  also  down  by  6  db  compared  to  3  db  for  the  noise  pre¬ 
ceding  the  signal.  Thus  a  loss  of  at  least  3  db  is  ex¬ 
pected  if  the  signal  model  at  APOK  is  taken  as  a  normal  up- 
gcing  pulse  and  echo.  For  the  imaged  array,  the  frequency 
of  F-K  spectral  peaks  appear  to  be  more  consistent  between 
the  signal  and  its  coda.  The  imaging  techniques  appears  to 
eliminate  conversions  other  than  a  simple  echo  at  the  surface. 
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Figure  1.  Simulated  signal  using  acoustic  log 
propagation  velocities  measured  at  APO1' 
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Figure  2.  Simulated  noise 
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Figure  3.  Ambient  noiBe 
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Figure  4.  Noise  sample  before  the  Aleutian  Earthquake 
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Figure  5.  Main  pulse  of  the  Aleutian  Earthquake 
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Figure  6.  Coda  following  the  main  pulse  of  the 
Aleutian  Earthquake 
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Figure  7.  Mirror  imaged  ambient  noise  sample 


VFKSPTRM  AMBIBIT  NOISE,  STATION  APOK 


SF|SHrUM»M  *0,  •  6  3>3 

NO,  or  CHANNEL  •  1» 

(MIRROR  IMAGED) 

01HHLINU  A. Ip  ■  an . U. 

ST  ART  I NQ  POIBT  • 

1  TOTAL  POINTS  .  »1» 

Ink  NUMBER  U*  SMOOtH|BC 

TlMft  •  * 

CHANNEL  It! 

SCALE  FACTOR 

DE2TH 

D  8 

Dwl 

1.06 

•2.91* 

0  -  •» 

UW2 

1.06 

•2.290 

A  -  V 

DuJ 

1.00 

-1.99* 

12  •  19 

Uw4 

1.00 

•1.666 

IB  -  21 

UuB 

1.00 

•  •015 

24  -  27 

Uu5 

1.C0 

.015 

Dw4 

1.00 

1.660 

UW3 

1.60 

1.990 

Owf 

liBB 

2.290 

Dwl 

1.00 

2.910 

•  A  •  i  N  II  •  •  I  * 

l.lll*  —,«.«» 

-t-  I  •.nii<  i  •  I  •*»••• 

i  ...  i  .  > 

i  t  •  i  .  i  ••••» 

i  i  .  i  ••  J  •••• 

i  i  •  ••«» 

,«**•••*  as ••••• *41 ••••••••• •*"*"•* •"** 

1 1 . lit  HI  lllll 

i  ••••■tM  ••••••••  |Ham  ■  i  ii— ■  > 

•s 

•miHi  I  i  itt>a 

in~r 

nr  -  660B0B6800 Xil/  • 1 11 

1  SMM  !  Ml  Ml 


Tu«  MAI*  ViCU«  • 


l.4mi  •• 

.Mil  .Ml*  .6*60 

I.....  ,  . . * - 

I  .  i  »  - 

I  I  I  ...  I 


•  •••  . 


r 

b 

E 

0  1.50 

tt 

E 

N 


•  .  lyiHM  M  • 

m  IIMMlIB  I 

. .  •  I  I 

. . I. . . 

I  l.l — •••IT -  ’  1 

•  •  •  •  *  .  «  2 

:  i.».  I 

t  ....  ••• 


»t..«  i  ...  mi  »••  i 

■  wk  mm i  •  2U22  i  it  i 

ttuutr  j/*n 

i  i  mS^ml  219m?  [»<l 

^rvtniiniir—»**?t*"V»5 

•  I  **"  ■  BBS  I  Ml 

Ml  BBB  I  B88B  1-  BBBB  I 


|  BBBBBB  6B  V  BB  B006B0 

. .  t  ••  B6BBB  m  I  .... 

♦ . .  •  ••  —  . . *”--.— 

|  BBB  i.t  B  •••  IB* 

|  ...  B  »*  BBIBB  ..  B  ... 

I .  B  •  BBBIBBB  :  B  .... 

.  I  ..  BBIBB  ..  I  .... 

*•.1.  »*i  I  00060  ...  I  •• 


mnhiuiwffi«MM,t**iltt*tM*  »i 

JGTYS.  I  II 

h*lJ*77??2922772?97?7)  I  0066666 I 

\  ^IbBBBBBBBB 

BBB  IBB  1  BBBBB  I  •! 

ABB  •  SB  BBMBftB  •••»•.  *i..»l 

I  *80860600  J 

...••I  I  . .  *1 


I - f.i.r*tt»»l 

I  I  .  I 

I  I  •  «  •  *  .4# • • • •  • 


*.•**••  it.  I 

......  ..  I 

.......  .  BBB 

■ . 000001 


.......  I  *'i.  .a.  I  .......  i 

•  •••••  ,,—♦••••••—  i .  j.-jir1”* . . 

irmtMi  I  I r»  I  •  ••••  I  •**  1 

I  ......  I  I  ...  !•.*/%••  I  •••  • 

|  I  I  ... ..«••«•/ I*.  I  ••  • 

I  I  I  . .  I  ••  • 

i.  I  or  i  •«•••»  ••••*•  i  .  . 


I  I  * .  /  I  •  •  (  (I  •  BBB  • 

|  |  .......... f  I  |  •  .  BBB  . 

I  «.|  |  •*•••*  ...V*.  L  •  OM  . 

* . .  3giia<«MN*i - 

I  n«  — rn«  •  •  ••/If* 

I  .......  |  .«••  I  •  *••!  •  ••!  I  f* 

I  f  ••  I  ••••  Ml  •  BilL^* 

I  I  •••*!•••  •«  ▼  »•  •  MN  • 


lllll  .............  I  .a  .... 

Qg;.::::.!..::;. — : — 

rl.ll  .1.  I  . .  •  ii  I  i*/ 

IBBBB  •  M.  |  .  .  ||A*  .  >1  ...  I 

BBB  .  i.  |  i.il .  I 

BBB  .  ..  |  ..1  .  t 

mb  ,  ,  i  . . »..  i 


I  ...V.i.fl 


BBBIBBB 

B  I  B 
B— ♦— 0- 


««  I  Ml 
•  I  *  i 

BBBB 

BBBBBB* 

nOjn 

--68*066- 


-f—  I  "nt«i»l  •••  I  ••  ••  B  ••  ••  I  *•• 

l  I  .*•••*••,•.••••••»••••  I  •  •  BBB  .  .  « . 

•  •••••4*  1  •!♦••••  BMB  •  BBBBB  .  BBBB 

. I  ••  ••Hl»#  •  BBBIBBB  .  BOMBBS 

I .  f  I  ••<!>•  .bib.  «*Q>bo 

. * . . . B0BBBB--.-B—  ♦— B-.—  BBBlTBB  — ' 

•  -  Itiui - -0*6*0  B6B6BB6  BBB  I  666  2**|*2*  BBB  I  BBS 

f  BB.  t  B  |*2?2?2*222*  22*12**2*  BBBBB  22222(22* 

ft  B  6BBBBB6B-  2*»*  **22  22  *  22  BB  |  BB  22  I  *2 

222  8  BOBBBBBBB  ***£  |  22*  2  600  266  |  BB2  6*6  * 

i  }2  BBBB  I  BBBBB  22f  *0  |  22  BBBBB  l/fl\  l  BBBB* 

•  — 2«BB8B0B«*B*e.2*-*— ABB0BBBB6  — 2— BB*BB0B--/*O***V--BB*B*BB  — : 

• — J-  66600*406*6-  2**( - BABB 44 B  2  2  |4|B|B  B  VllUt  BBBBB*  * 

222  6  I  6  ***  I  92  *22  B  JBB^T^BB*  B  **« 

•  BB. . .  22*22*2*2  B  ***(  2*2  66(66  22*  1*22 

BB  ,  |  ..  (BBbBaBBBBBBB  •  BB  I  BBB  22  I  22  BBB  I  BB  . 

|  „»•  I  h;  I  B  .6  2  6.  B  i. 


I  I 

. . . 

D6BBBBB  BBBBB 
22*22222*22! 
*2*2  2* 

*2*  |  *2 


t  t  I 

.  I  t 

I  . I 

I  . I 

I  . I 

-♦ . . 

BB  ,~;i..l  ■  I 

2f  B  I  .BB  I 

2*22  BBBBBBBB  *  21 

****  BBBbBBBBB  B  22| 


*2  |  BB  2*2  BBBBBtBBBB  2*  I 

2— BBBBBBBBB— ♦-22.000000B0000.2— ♦ 
*  6666*60  12*2  00006660000  *  | 
*21  222  B  I  B  2*1 

8  222*2222*  |  BB  . BB  | 

BBBB6B6B686BI  I  .  M 

• .  I  ...  I  ••  U 


ft  |  '.PMWwmu'  <■.  *•***•  «**•*  •••••«*  I  •  tl  *»• 

*2  |  6*66666666  2*2  606666  1*080  BBBBBB  262  6666666665  I  2*  BB| 

—  »  A.AAAAAAA*  ****  A.AAAA  ••**••  **•■»  AAA.......  .  <*-  ••• 


•V  'IT 

B*  B  | 


Figure  8.  Mirror  imaged  noise  before  Aleutian  Earthquake 
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Figure  9.  Mirror  imaged  main  pu  1  a'e  of  Aleutian  Earthquake 
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Figure  10.  Mirror  imaged  coda  of  Aleutian  Earthauake 
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1 11.  ABSTRACT 

The  seismic  date  analyzed  include  an  Aleutian  earthquake,  several 
samples  of  ambient  noise  along  with  the  noise  occuring  before  the  event,  and 
a  25  second  sample  of  coda  following  the  signal.  The  signal  coda  sampled 
had  nearly  uniform  power.  A  preliminary  Interpretation  of  the  noise  and 
signals  strongly  suggests  the  importance  of  energy  conversions  possibly 
due  to  the  complex  geology  at  APOK.  The  geology  is  known  to  be  very  complex 
and  is  characterized  to  first  order  a  thick  wedge  of  low  velocity  strata 
'abedded  in  material  of  higher  propagation  velocity.  The  evidence  of  con¬ 
version  is  based  on  the  asymmetry  of  the  T-K  spectrum  of  ambient  noise  with 
re', pact  to  positive  and  negative  wave  number.  The  resultr  indicate  a 
predoi i nance  of  obliquely  incident  up-going  waves  at  all  frequencies.  Further, 
the  signals  and  coda  similarly  displayed  anomalous  down-going  pulses  atten¬ 
uated  by  approximately  6  db  with  respect  to  the  up-going  pulse.  Further, 
the  down-going  pu.Vse  contains  lower  frequencies  than  the  up-going  pulse  and 
has  a  much  lower  apparent  vertical  velocity.  The  results  obtained  by  mirror 
imaging  the  vertical  array  indicated  that  a  loss  of  approximately  3  db  would 
be  encountered  in  using  a  conventional  signal  model  at  this  site. 
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